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Abstract 

We calculated bulk viscosity due to non-equilibrium weak processes in superfluid nucleon-hyperon 
matter of neutron stars. For that, the dissipative relativistic hydrodynamics, formulated in paper [lj] 
for superfluid mixtures, was extended to the case when both nucleons and hyperons are superfluid. 
It was demonstrated that in the most general case (when neutrons, protons, A, and S~ hyperons 
are superfluid), non-equilibrium weak processes generate sixteen bulk viscosity coefficients, with 
only three of them being independent. In addition, we corrected an inaccuracy in a widely used 
formula for the bulk viscosity of non-super fluid nucleon-hyperon matter. 

PACS numbers: 97.60.Jd, 26.60.+C, 47.37.+q, 47.75. +f 
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I. INTRODUCTION 



It is well known (see, e.g., , that neutron stars can be unstable with respect to 

the emission of gravitational waves. The matter in a pulsating neutron star is not (even lo- 
cally) in chemical equilibrium. The relaxation towards chemical equilibrium is accompanied 
by the dissipation of pulsation energy. This process is one of the most important dissipative 
processes, suppressing the growth of gravitational-wave instability. It can be described by 
the introduction of an effective bulk viscosity in the hydrodynamic equations (see, e.g., 6j). 

The bulk viscosity due to non-equilibrium weak processes in neutron stars was calculated 



by a number of authors (for a review see, e.g., |7j). It strongly depends on the composition 
of stellar matter. In this paper we consider the matter of inner layers of neutron stars, 
composed of electrons, muons, neutrons, protons, A, and S~ hyperons (nucleon-hyperon 
matter). A calculation of the bulk viscosity for nucleon-hyperon matter is complicated by 
the fact that baryons in such matter can be superffuid Q, [of, \ 
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The bulk viscosity of superfluid matter was calculated in a number of papers (see, e.g., 
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17|). In these papers only one bulk viscosity coefficient was studied, 



analogous to that in non-superffuid hydrodynamics. The effects of superfluidity were taken 
into account only in calculating the reaction rates. However, it is known that there are 
several bulk viscosity coefficients in hydrodynamics of superfluid liquid [l|, E3, Ji], 19]- 

In a recent paper pj four bulk viscosity coefficients were calculated for the matter com- 
posed of superfluid neutrons, superfluid protons and electrons. It was shown that taking 
into account three additional bulk viscosity coefficients results in a significant decrease of 
characteristic damping times of sound modes (approximately, by a factor of three). 

In this paper we extend the results of Ref. l( to the case of nucleon-hyperon matter. 
In particular, we demonstrate how the dissipative hydrodynamics [II] should be modified 
to describe a possible presence of superfluid hyperons. Next, we show that in the most 
general case, when baryons of any species are superfluid, non-equilibrium processes of mutual 
transformations of particles generate sixteen bulk viscosity coefficients, only three of them 
being independent. In addition, we correct an inaccuracy in the expression for the bulk 
viscosity of non-superfluid nucleon-hyperon matter made in Ref. [15| and spreaded widely 
in the literature (see, e.g., l2ll . |22| . |23| . 124 l25ip . We calculate the bulk v scosity correctly 



and compare our results with those of Ref. 
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The paper is organized as follows. In section II we calculate the bulk viscosity of non- 
superfluid nucleon-hyperon matter. In section III we calculate and analyze all sixteen bulk 
viscosity coefficients describing dissipation in superfluid nucleon-hyperon mixtures; the re- 
lations between these coefficients are also discussed. Section IV presents a summary. 



II. BULK VISCOSITY OF NON-SUPERFLUID HYPERON MATTER 



In this section we derive an expression for the bulk viscosity due to non-equilibrium 
processes of particle transformations in a dense non-superffuid matter composed of electrons 
(e), muons (//), neutrons (n), protons (p), and hyperons (A and S~ hyperons). Here and 
below the variation 5 A of some physical quantity A will be defined as the difference A — Aq, 
where Aq is the value of A in thermodynamic equilibrium (when matter is unperturbed). 

The most effective weak processes in nucleon-hyperon matter are the following non- 

nrr 



leptonic reactions 

mm 
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27] 



n + n <-> p + £~, (I) 

n + p <-> p + A, (2) 

n + n <-> n + A, (3) 

n + A <-> A + A. (4) 

The full thermodynamic equilibrium implies the equilibrium with respect to these reactions, 

2/X n o - /ipO - /"SO = 0, (5) 
A*nO - Hao = 0. (6) 

Here //jo are the chemical potentials of particle species % = n,p, A, S taken at equilibrium 
(correspondingly, fa are the chemical potentials in the perturbed matter). Notice that, the 
equilibrium conditions for reactions (TSJ), (j3J), and (jlj) coincide. 

Leptonic reactions (e.g, direct and modified Urea processes with electrons or muons) are 
much slower in comparison to the reactions ([!])- (jlj). For 'typical' perturbation frequencies 
(e.g., 10 3 — 10 4 s _1 for radial modes or for r-modes of rapidly rotating neutron stars) the lep- 
tonic reactions cannot influence substantially the chemical composition of perturbed matter. 
Hence, the main contribution to the bulk viscosity comes from the non-leptonic reactions 
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CO)-®- In addition to the processes described above, there is a fast non-leptonic reaction 
due to the strong interaction of baryons 



n + A <-> p + S . (7) 

In accordance with Ref. [15] we assume that the perturbed matter is always in equilibrium 
with respect to this reaction, 

<Wst = fi n + A*A - - A*s = 0. (8) 

Let us obtain the expression for the bulk viscosity of non-superffuid nucleon-hyperon matter. 
For that, we consider a pulsating nucleon-hyperon matter, slightly perturbed from an equi- 
librium state (so that one can use the linear perturbation theory). If the reactions (jTj) — (J^J) 
are forbidden, then pulsations are reversible and there is no energy dissipation (notice that, 
the reaction [7] is open). We denote the pressure in this case by P cq . The presence of the 
reactions (jXJ) — (SI) in the pulsating matter leads to a difference between the real pressure P 
and P eq . We define the bulk viscosity £ by the formula 

P-P cq = -£div(«), (9) 

where u is the hydrodynamic velocity of pulsations. Notice that, this definition differs from 
the usually accepted one (see, e.g., Ref. [l5]). Usually, instead of P cq in formula ([9]) it 
is common to substitute the pressure which would be established in the pulsating matter 
assuming that there is an equilibrium with respect to all the reactions (i.e. the reactions are 
very fast). Both these approaches are possible. 

Generally, the pressure P and the other thermodynamic quantities depend on six pa- 
rameters, for example, the number densities rij, where j = n,p, A,S,e, /i (one can neglect 
;he dependence on temperature in strongly degenerate neutron-star matter, see, e.g., Refs. 



II, 128|, I29j). However, these parameters are not all independent because in the nucleon- 
hyperon matter two conditions should be satisfied: the equilibrium condition (JHJ) with respect 
to the reaction (j7j) and the condition of quasineutrality, 

n p = n e + n^ + riY,. (10) 

Taking into account that the reactions (flj)-(pi| and (JTj) conserve the number of leptons and 
the leptonic processes are neglected, we obtain that the relative number densities of leptons 



x e = n e /n b and x^ = n^/n b (n^ = n n + n p + n\ + ns is the baryon number density) remain 
constant during the pulsations, 

5x e = 5x^ = 0. (11) 

This result is valid only for non-superfluid matter and follows from the continuity equations 
for baryons, electrons, and muons, 

^ + dw(n b u) = 0, (12) 
^ + divKn) = 0, (13) 
^ + div(n M «) = 0. (14) 

If baryons of any species n, p, A and/or £ are superffuid, then the continuity equation for 
baryons (fl2"j) should be modified (see Section 3) and Eq. (ITT1) does not hold. 

In view of Eqs. (jBj) and ({TO]) , pressure is a function of only four independent vari- 
ables, say, n b ,n H ,x e , and x M (n^ = ba + is the hyperon number density). Expanding 
P(n b , x ei %fj.) in the Taylor series near the equilibrium state, one obtains for the variation 
of pressure SP, 

dP(n b ,n H ,x e ,x^) dP(n b , n H , x e , , . 

5P = — f)n fe H —on H) (15) 

on b ariii 

where we used Eq. (II II) . The variations Sn b and Snn can be found from the continuity 
equations for baryons ffl~2]) and hyperons, 

d5n H 



dt 



+ div(n H u) = ATi + Ar 2 + Ar 3 + Ar 4 . (16) 



Here Ari, Ar2, Ar3, and AT 4 are the net numbers of hyperons generated in a unit volume 
per unit time in reactions ([Tj), d2J), and (j4]), respectively. If deviation from the equilibrium 
state is small, the sources AI^ (1=1,. . .,4) can be expressed as (see, e.g, 
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IM) 



Ar, = XJhi, (17) 

where A; are the 'reaction rates', some functions of number densities and temperature; 
8 fa = 2/i„ — n P — 5/J.2 = S Hz = <5 Hi = A*n — A 4 a are the chemical potential disbalances 
for the reactions ([I]), (j2J), (EI), and (jlj), respectively. Taking into account the equilibrium 
condition (jHJ) for the fast reaction (JTj), one has: 5 Hi = &Hi = ^3 = ^A*4 = ^A 4 - 
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Notice that, there is no source in the equation (ITS]) owing to the fast reaction (JTj), because 
this reaction does not change the number of hyperons. The choice of another variable instead 
of nu (for example, the neutron number density n n ) would make it necessary to take into 
account the source due to the reaction (j7| 



the subsequent papers [20, 



21 
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23, m, 



n the paper of Lindblom and Owen [15J (and in 
251 ]) the number density of neutrons was chosen as 
such variable, but the source of neutrons owing to the fast reaction (JTj) was neglected. This 
leads to an error in the expression for the bulk viscosity. 

In fact, one could think that the source of neutrons due to the fast reaction (J7|) equals 
zero, because the matter, as we mentioned before, is in equilibrium with respect to this re- 
action. However, this is not quite true, because even small (negligible in all other situations) 
deviation from the equilibrium <5/if ast = fi n + A*a ~ A* P ~ A*s multiplied by the large reaction 
rate Afast results in a finite (non-zero) source Arf ast = Af as t^fast- This fact was emphasized 
by Jones [261 ] . 

Now let us assume that the perturbation of matter is periodic, so that all the thermo- 
dynamic quantities oscillate near their equilibrium values with the frequency uj. Then one 



finds from the continuity equations ( |T2l) and ([IE 

8n b = - 



— div(tt) 

iuj 



8n H = — — [nffdiv(u) — XSfj], 



(18) 
(19) 



where A = A1 + A2 + A3 + A4. As in the case of the pressure, the chemical potential disbalance 
8\i in Eq. ( [191 is a function of n b , %, x e , and x M . In analogy with Eq. ( TTBT) for 5P, it can 
be expanded near the equilibrium state and written as 



d5fi(n b ,n H ,x e ,x ll ) x d8/i(n b , n H , x e , x^) 
on b + 

on b oriii 



5n H . 



(20) 



Here we take into account that 5fi = in the full equilibrium and that the relative number 
densities x e and x M do not change in the course of pulsations (see Eq. (TTTh ) . 

Now, solving the system of equations (jT8"]) - (|2T)]) and substituting the expressions for Sn b 
and Sn H into Eq. ([TBI) , we derive 

dP{n b ,x H ,x e ,x^) A d8fi(n b , x H , x e , x M ) 



P-P, 



cq 



div(tt 



x 



8xh oj 
iX d8fi(n b ,XH,Xe,XiJ.) 



dn h 



um b 



dx 



+ 1 



-1 



(21) 
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Here the independent variables are n b , xh = nn/n b , x e , and x^. It is easy to express the 
bulk viscosity £ from this equation. We are mainly interested in the real part of £ because 
it is Re(£) that is responsible for the energy dissipation (see, e.g., [121] ). In this sense it is 
probably more appropriate to define Re(£) as the 'real' bulk viscosity. It equals 



n\ dP(n b ,x H ,x e ,x^) d6fi(n b , x H , x e , x^) 

Ree = "T dx~ H ^ 



86lJ,(n b ,XH,Xe,Xn) 



dx 



H 



1 + U 2 T 2 



(22) 



where r = n b /X [d5fi(n b ,XH,x e ,x^) / 'dxn} 1 ■ 

For comparison, we present here the result of Lindblom and Owen |15J (notice that, these 
authors neglected the reactions [3] and H[ thus assuming that A3 = A4 = 0) 



Re £ L = - 



n b dP(n b , x n , x e , x^) d5fi(n bj x n ,x e ,x^ 



2Ai + A 2 



dx r 



dn h 



d5fi(n b ,x n ,x e ,x^ 



1 -2 



dx r 



(23) 



where r L = n b / (2Ai + A 2 ) [dSfi(n b , )/dx n 



X 1 ■ 



n n /n b . 



The bulk viscosity ([22]) depends on the reaction rates Ai, A 2 , A 3 , and A 4 . Some of them 



were calculated in a number of papers 
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271 ] . Different authors present different 



results for the rates; the discussion of advantages and disadvantages of their calculations can 
also be found in those papers. Unfortunately, to date, there are no strict calculations of the 
reaction rates. It is reasonable to think that all the rates are of the same order of magnitude. 

The dependence of the bulk viscosity on the baryon number density for the temperature 
T = 3 x 10 9 K and the oscillation frequency u = 10 4 s -1 is presented in Fig. 1. While 
calculating the bulk viscosity, we used the third equation of state of Glendenning 301 ] . 

The solid line illustrates our results for the bulk viscosity obtained from Eq. (!22j) . We 
employed the reaction rates from Ref. [ijj] and, following that paper, we put the rates of 
the reactions ([2D and (@J equal zero, A 3 = A 4 = 0. The dot-dashed line is the bulk viscosity 



calculated as described in Ref. 



151 ] (see also formula 123]) . We remind that in that paper 



the relative number density of neutrons x n was chosen as one of the independent variables. 
However, the source of neutrons due to the fast reaction ([7]) was erroneously neglected. 
As one can see, this mistake does not influence the results si gnifi cantly (typically, by ~ 
(10 - 30) %). By the long dashes we show the results of Ref. 14(. In Ref. [14j, only one 
hyperon reaction was taken into account, namely the reaction ([T]). As in Ref. [l5j, the source 
due to the fast reaction ([7]) was neglected and, in addition, the equilibrium condition (jSJ) 
with respect to this reaction was ignored. Moreover, the authors of Ref. [lij] used the non- 
relativistic approximation when calculating the rate of the reaction (00). This assumption is 
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FIG. 1: Bulk viscosity £ versus baryon number density at T = 3 x 10 9 K and w = 10 4 s 1 
for non-superfluid matter. Solid, long-dashed, and dot-dashed lines show our results, the results 



of Ref. 



1J and Ref. 



respectively. Vertical dashes indicate the thresholds for (from left to 



right): appearance of muons; direct Urea processes involving electrons and muons, respectively; 

y. The inset demonstrates the difference between 



appearance of A and S hyperons, respective 
our calculations and those of Refs. [14| and 



151 ] in more detail. 



not well justified for the baryon number densities in the range rib > (0.3 — 0.6) fm 3 because 
of the strong dependence of the reaction rates on n& (see, e.g., Ref. [151]). 



III. BULK VISCOSITY OF SUPERFLUID NUCLEON-HYPERON MATTER 

In this section, unless otherwise is stated, the subscripts i and k refer to baryons (i, k = 
n,p, A, £). The summation is assumed over repeated baryon indices i and k. The subscript / 
refers to leptons (/ = e, /z); the subscript j runs over all particle species (j = n,p, A, E, e, /i). 



A. The relativistic hydrodynamics of superfluid nucleon-hyperon mixture 

In Ref. [l| the dissipative relativistic hydrodynamics of superfluid mixtures was formu- 
lated for npe matter. Here we extend this hydrodynamics to the case of superfluid nucleon- 
hyperon matter composed of superfluid protons, neutrons, A and S~ hyperons, as well as 
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normal electrons and muons. 

The general formulae (Eqs. 26-34 of Ref. [lj) describing the relativistic hydrodynamics 
of superffuid mixture remain valid with the notion that now the subscripts % and k refer not 
only to superffuid nucleons (i, k = n,p) but also to superffuid hyperons (i, k = n,p, A, E). 
For instance, the continuity equations for particle species j are written as 



with 



Energy-momentum conservation law has the form 



(24) 



(25) 



d^ v = 0, 



(26) 



where 



T iw = (p _|_ e ) u n u » + Pr] 



[IV 



+Yi 



ik 



V 



(27) 



and t^ v is the dissipative correction to the energy-momentum tensor which will be specified 
below. It satisfies the constraint 

u^T^ = 0. (28) 
The hydrodynamic equations must be supplemented by the second law of thermodynamics 



Y 



ik 



de = T dS + fii dn, + ji e dn e + /x M dn^ + — d 



(29) 



Using the quasineutrality condition ( TTDT) and the condition of equilibrium (jSJ) with respect 
to the fast reaction ([7]), Eq. (1291) can be rewritten as 



de = TdS+fin dn b -5iidn H -(lJ>n-fJ>p-fJ'e) dn e -(fi n -fi p -fi^) dn^ + -^-d 



(30) 



where we remind the notation 8[i = fi n — In full thermodynamic equilibrium the third, 
fourth, and the fifth terms are zero because of Eqs. fl5J) and ((6]), and of the beta-equilibrium 
conditions, \x n = \x v + pi e and ji n — ji p + fifi (see, e.g., Refs. |l2j, If3l|). 

In Eqs. ( 124")) - (130]) Y ik is a 4 x 4 symmetric matrix which is related in the non-relativistic 
limit to the entrainment matrix p ik by the equality Yik = Pik/(mim k ), where m» 



is the mass of a free baryon of species i (the matrix is a natural genera. 
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32 



ization of the 
33j). To the 



superfluid density to the case of superfluid mixtures, see, e.g., Refs 
best of our knowledge, the matrix has not been calculated for a nucleon-hyperon matter. 
Furthermore, e is the energy density; fij is the relativistic chemical potential of particle 
species j; and S is the entropy density. The pressure P in Eq. (|2"T|) is defined in the same 
way as for ordinary (non-superfluid) matter l|, |34( , 



-e + iiiUi + fi e n e + + TS. 



(31) 



Next, rf" = diag(— 1, +1, +1, +1) in Eq. (l27j) is the special relativistic metric; is the 
four- velocity of normal (non-superfluid) liquid component normalized so that u^u^ = — 1 
(we assume that all non-superfluid components move with the same velocity u^). The four- 
vectors wf-s satisfy the condition 







(32) 



and describe motion of superfluid components. To take into account the potentiality of 
superfluid motion, a four-vector u>^ should be expressed through some scalar functions fa 
and written as (see Ref. [l[) 



wfc = d»fa - qi A» - fa + k,)^. 



(33) 



Here the scalar fa is related to the wave function phase of the Cooper-pair condensate $j 
by the equality V0j = ^V$j/2; is the four-potential of the electromagnetic field; is 
the electric charge of particle species i; Xi is a small dissipative correction to be determined 
below. 

Note that one can avoid the introduction of new functions fa in the hydrodynamics of 
superfluid mixtures if one formulates the potentiality condition fl33l) in the equivalent way 



(34) 



Below we will use the latter formulation because it is more suitable for our purpose. In this 
approach, four-vectors are treated as independent hydrodynamic variables. 

The hydrodynamics discussed above would be incomplete without an indication what 
we mean by a comoving frame, that is the frame where we measure (and define) all the 
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thermodynamic quantities. As was demonstrated in Ref. [lj], the condition (|32|) dictates 
that the comoving is the frame where the four-velocity equals = (1,0,0,0). In this 
frame, the basic thermodynamic quantities s, rij, and (or V0j) are defined by (see Eqs. 

(ESD, (EZD, (EHD, and (P5) 



h = n i> 

rpOO 



Y ik W( k) = Y ik V(J) k} 



e. 



(35) 
(36) 
(37) 



All other thermodynamic quantities can be considered as their functions or, equivalently, 
the functions of e, rij, and w'tiW^p- 

In analogy with Ref. [lj], from Eqs. ( 1241) . ( |26i) . and (l29i) one can derive the entropy 
generation equation, which defines the dissipative corrections and Xj, 



Here is the entropy current density, 



(38) 



(39) 



satisfying the natural constraint u M S* M = — S. The last two terms in Eq. (1381) are small and 
can be omitted in the majority of applications (for more details, see the discussion in Ref. 

Using the requirement that the entropy does not decrease, one can obtain for the dissi- 
pative corrections (neglecting the last two terms in the right-hand side of Eq. 



-K 



X, 



/' 



(H™ u u + u") (<9 7 T + Tu s d 5 uj) 
iuH^d, \Y lk wl^\ -i 2 H^d lU \ 

- Up <V* M . 



-in d» 



Y ik w^ k) 
Y ik wfo 



(40) 
(41) 
(42) 
(43) 
(44) 
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Here k and r] are the thermal conductivity and shear viscosity coefficients, respectively; 



H* w = rj^ u + u^u u is the projection matrix; ( G- G>- £w- £•">;• an d £7* are twenty five 
bulk viscosity coefficients (i = n,p, A, E). We would like to emphasize that the dissipative 
corrections (|40j) - (|44p are incomplete in a sense that they (i) do not include various terms 
related to particle diffusion; (ii) neglect (typically) small terms, explicitly depending on w'L. 
For example, we neglect the terms of the form w^d^T and u^w^d^u 1 in the expressions for 
Xk and r^, respectively. The similar approximation is very well known in the non-relativistic 
theory of superfluids (see, e.g., the textbook of Landau and Lifshitz Q, §140orRef. Q). An 
inclusion of all these dissipative terms would lead to a number of kinetic coefficients much 
larger than 27. Since in this paper we are mainly interested in the bulk viscosity coefficients, 
we restrict ourselves to a simplified form (|40p - (|44p of dissipative corrections. 

The number of bulk viscosity coefficients can be reduced. Notice that, from the quasineu- 
trality condition (TTOT) and the continuity equations (T24|) . it follows that 



Y 



p kW (k) 



(45) 



A similar condition for superfluid npe matter was derived in Ref. I] (see Eq. 41 of that 
reference). Owing to the condition (145 p . there is no need to introduce the bulk viscosity 
coefficients both for protons and for S~ hyperons. For example, it is sufficient to introduce 
the quantity Gsp = £i p + £xs instead of £i p and Because of the same reason, we are 
not interested in the quantities x^ and x p taken separately. Instead, we will introduce the 
sum xs p = xs + x p (notice that, xs p is the quantity that appears in the entropy generation 
equation 1381) . As a result, the corrections t^ u and x g take the form 



-k (H^ u v + H ury u") (<9 7 T + Tu s d s u 7 ) 
7] H" 5 \ d sUl + d^us - ^ v-ys d £ u £ ^j 



X A 



Y 



7 

7" i 



1 kW (k) 



6 d jU 



Y qk W( k) 



Y 



- Ga 9^, 



(46) 
(47) 
(48) 
(49) 



Here and below the subscript q takes on the values n, A, and Sp. The expression 



Y 



T,pkW {k) 



in Eqs. p6|)-(l49j) means d„ 



Yvkwfa 



or, equivalently, du 



Y, 
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As follows from the above equations, we have actually sixteen (rather than twenty five) 
bulk viscosity coefficients which can contribute to the dissipation of mechanical energy in 
superfluid nucleon-hyperon matter. Using the Onsager symmetry principle, we obtain 



Thus, generally, only ten of them are independent. 

B. Calculation of the bulk viscosity coefficients for superfluid nucleon-hyperon 
matter 

Let us calculate these phenomenological coefficients assuming they are due to non- 
equilibrium processes (|T])-(jIJ). As in Sec. II, we assume that the matter is slightly perturbed 
out of equilibrium state and pulsates with the frequency u. Since the deviation from the 
equilibrium is small the hydrodynamic equations can be linearized. Because of the same 
reason, the dependence of various thermodynamic quantities (e.g., the pressure P) on the 
scalars w^w^ti can be neglected. We assume that the four-vectors w?* characterizing the 
superfluid flow of particle species i = n,p,A, or £ equal zero in the unperturbed matter. 
Thus, by perturbing the system, we produce some small so that the scalars w'Lw^^ 
will be of the second-order smallness and can be omitted in the linear approximation. 

We start from the non-dissipative relativistic hydrodynamics of superfluid nucleon- 
hyperon mixture. In this case the energy- momentum tensor is given by Eq. (1271) . where 
the dissipative correction should be set zero. Similarly, the expressions for are given 
by Eq. ( 1341) with Xj = 0. The non-equilibrium processes (jTJ) — (SI) lead to the appearance of the 
sources in the right-hand side of the continuity equations (|24p which, as we will demonstrate 
below, generate the 'effective' dissipative corrections t^ v and x ? . 

To calculate £\ q and £2 it is convenient to expand the energy-momentum tensor of the 
pulsating matter (1271) (with t^ u = 0) in the comoving frame (where = (1,0,0,0)) near 
the equilibrium, as it was done in Ref. [lj], 



<^3A — ^5n, ^3Sp — ^6n, ^4n — £ln, ^5Sp — ^6A; &A — £lA, ^4Sp — SlEp- 



(50) 



T 



•00 



e + 5e 



r 



T 



in 



(fc)> 



T 



(P + SP) 5; 



(51) 
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Here we restrict ourselves to the linear perturbation terms. The spatial indices n and m vary 
over 1,2, and 3; Eq, fao, and Pq are the corresponding thermodynamic quantities calculated 
at equilibrium (in the absence of pulsations). 

Now our aim is to extract from the tensor floTl) various terms which are generated by the 
non-equilibrium processes and contribute to dissipation. Then, we will write these 

terms in the form of a separate dissipative tensor r^ lk . 

As follows from Eq. fl30l) . in the linear approximation 5e remains the same as in the 
absence of dissipation, 5e = [i n 5rib (this is because the dissipative processes (p])-® conserve 
the number of baryons, hence Srib is independent of the reaction rates A;, / =1,2,3, or 4). 
Thus, the component r™ lk of the tensor r^ lk is zero. Similarly, r^ k = r^ k = 0. On the 
contrary, the variation SP of the pressure contains a dissipative part (we denote it by 5-Pdis)- 



According to Refs. 
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15| and Sec. 2, it is given by 

5P dis = Re(5P), (52) 



so that the tensor r^n, can be presented in the form (in the comoving frame) 



r 00 - n 

r bulk _ u > 

T-Om _ mO _ n 
'bulk — 'bulk — u ' 



CTk = ^dis <W (53) 



Let us determine 5-Pdis- For that purpose, we present the pressure P as a function of n b , 
fin, ""-En = nj] + n n , and y = Xe/x^. All other number densities can be expressed through 
rib, nn, n-£ n , and y with the help of Eqs. (JH1) and (|T0l) . Notice that we choose ns„ and y 



instead of the variables x e and x^ of Sec. 2. The variables x e and x M are less convenient here 
because Eq. (fTTT) does not hold in the case of superffuid matter. Expanding the pressure 
P(rib,nH,nzn,y) near the equilibrium state, we write 

dP dP dP dP 

SP dis = -w- Re(Sn b ) + ^ — Re(5n H ) + Re(5n Sn ) + — Re(^). (54) 

on b on H OriY.n oy 

It is straightforward to show that 

by = 0, (55) 

as a consequence of the continuity equations ffTBl and f|T4l) for electrons and muons. The 
variations of other variables, nb, riu, and ns n , can be found from corresponding continuity 



14 



equations. Using Eq. fl2^|) and the fact, that the variations depend on time t as exp(iu>t), 
we obtain in the comoving frame 



iuj 5n b + div(J b ) = 0, 
iuj 5riH + div(J#) = Ar, 
iu 5n Sn + div(J Sn ) = -Ar. 



(56) 
(57) 
(58) 



Here Ar = ATi + Ar 2 + Ar 3 + Ar 4 = A 5li, and 

Fit 

5n(n b} n H} nx n 



85ll 85ll 85li 

5n b + - — 5n H + 7; m Sn . 



8n 



8n 



8n 



(59) 



In Eqs. (ESD-dSBD Jb = n b u + Y,t Y ikW{k)] Jh = n H u + Y^ k w^ k) + Y Ak W( k y, J^ n = n^ n u + 
YskW(k) + Y n kW(k)', u and w^) are the spatial components of the four- velocity and four- 



vector w^y respectively. The solution to the above system of equations gives 



Re(5n b ) = 0. 
Re(<5n Sn ) = - 
Re(5n H ) = k 



Re(5n.ff) 
85 li 



dm 



div(J 6 ) + |^div(J H ) + i^div(J Sn ) 



8n 



H 



8n 



En 



(60) 
(61) 
(62) 



where we use the notations 1/k 



\(d5(j,/driH — d5n/dn^ n y (1 + tu 2 r 2> 



1/A (85Lt/8n H - dSfi/dnvn)- 1 . Now, using Eq. flUD for 5P dis and Eqs. p?j)-(|62l, we can 
find the dissipative tensor r^ lk in the comoving frame (see Eq. (153]) ). Then it can be easily 
rewritten in an arbitrary frame. The result is 



t; 



bulk 



H^k 



8P 



8n 



H 



8P 

dnj; n 



85/i 85/i ...... . 

n & h n H h n S n 7; <? 7 ii 7 

c>n b dn H 8riY,n J 



85 n \ 



an b 8n H J \ 



(63) 



A comparison of r^ lk with the phenomenological dissipative tensor r M1/ (see Eq. ( 1461) ) al- 
lows us to identify the bulk viscosity coefficients £i n , £ia, £ie p , and £ 2 > generated by non- 
equilibrium processes ([I])-® 
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£ ( 9P - ( d5 ^ + (64) 



' <9P 


dP 


<9n ff 


dn-zn 


' ap 


dP 




dnxn 


' ap 


dP 




dnxn 


dP 


dP 


dn H 


dnxn 



( dP dP \ ( db\x 88 /I 85 fi \ 

6 = -« q q 716-5 1" n HT> h TlSn^ 

Van H an Sn / V on b 8n H on^nj 



, , ( 67 ) 



It can be shown, that the expression for £2 formally coincides with Eq. (1221) for the bulk 
viscosity of normal matter (however, these formulae give different numerical results because 
the reaction rates A/ (I — 1,. . .,4) differ for superfluid and normal matter). To prove this, 
we can rewrite Eq. (|67|) using n b , xh, x e , and x M as independent variables instead of n b , n^, 
ris n , and y. The following equalities will be helpful (/ is an arbitrary function) 

1 df{n h ,x H ,x e ,x,j) df(n b ,n H ,nx n ,y) df(n b ,n H ,nx n ,y) 



n b dx H dn H <9n Sn 

df{n b ,x H ,x e ,x^) df(n b ,n H ,n^ n ,y) 

n b a = n b a 

on b on b 



(68) 



df(n b , n H , n Sn , y) df{n b ,n H ,n^ n ,y) . . 
+n H k N Sn 5 • (69) 

To calculate other bulk viscosity coefficients let us apply the same consideration to Eq. 
fl34"|) for w^y As a result, we obtain (in the comoving frame) the expression for the dissipative 
component Xj generated by nonequilibrium processes (pQ)-([Ij) 

Hi = 5ni dis = ^Re(5n b ) + ^-Re(5n H ) + -p^-Re(5n Sn ) . (70) 
on b on H on^n 

Here <5/Xjdi S is the dissipative term in the Taylor expansion of the chemical potential fii near 
the equilibrium state (similar to #Pdi S ); Re(Sn b ), Re(5n H ), and Re(5n Sn ) are taken from Eqs. 
(I60l) -( l62l) . In a fully covariant form, x g is given by Eqs. (l4"?l) - (l4"9~l) with the bulk viscosity 
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coefficients 



^3A 

^5A 
^6Sp 



-A; 



-A; 



—k 



-k 



dfJL n _ 8fi n 

8n H 

8fi n 



8n H 

8fi n 

8n H 
dfi n 
8n H 
dfiA 
8n H 
d{n P 



driT, n 

dfi n 
dn^n 

dfi n 

dfl n 

dn-zn 
dfiA 



dn-zn 
85 fi 



85 fx 85fx 
8n b 

d5fx 

dn b 8uh 

^85 fx 85 fx 85 fi \ 

8n b 8n H dn^n) 

85 fi 85 fx 

n b ^ hn/f7 hn En . 

8n b 8n H 8n^ n J 



85 fx \ 



n b 



8n^ n 
fix) d(fi p + /i S 



85 a 85 a 85fi\ 
8n b 8n H 0n Sn J 



8n H 
8fi A 8fi A 



-k 



-k 



-k 



-k 



8n En 

8fi\ 

dn^ n 

dfiA 



dn H 

dfiA 
dn H 

dfiA 

dn H dn^n 
d{fi p + /fg) _ 

dn H 
djfip + fig) _ 

8n H 
8{fi P + fix) 



+ 



dnx n 
85fi 
8n b 
85fi 
8n b 8n 



85 fi 85 fi 

n b -^ h n H - h n Sn - , . 

8n b 8n H dn En J 



85 fi \ 



85 fi 

dn^n 
85 fi 



H 



2 d5fi d5fx_ + d5fi \ 
dn b 8n H drix n J 



8(fip + fix) 
dn-zn 

djfip + fin) 
dn^n 

d(fi p + /is) 



dSfi 85 fi \ 

8n b dnxn J 

85 fi 85 fi 

8n b dnn 

2 d5fi 85fi_ 85fi \ 

8n b 8n H dn^nj ' 



(71) 
(72) 
(73) 
(74) 
(75) 
(76) 
(77) 
(78) 
(79) 
(80) 
(81) 
(82) 



8n H dn^n 

In Eqs. (1TTi) - (182l) we assumed that the thermodynamic quantities are functions of n b , Uh, 
TT-Sn; and y. One can easily check that all the six Onsager relations (j50p are satisfied. 
Moreover, it turns out that the bulk viscosities (!64|) - (l67|) and (iTT|) - (l8"2"|) obey a number of 
additional constraints (q = n, A, Hp) 



6 



f 2 

Sin 



1A 



£,3q + £5g, 6t£p = iin + ^4A, 

&£,3n, ClA = £lSp — S^Ep, 



(83) 
(84) 
(85) 



so that we have only three independent bulk viscosity coefficients, say £2>£i»n an d £ia- All 
other coefficients can be expressed through these three. The coefficients £2,6in> and £ia 
are compared in Fig. [2j Because the dimensions of the coefficients £ii {i — ^ ; A) and £2 
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N 1 



X 







n — i — i — i — i — r 



~i — i — i — r 



n A X £l A /^ 



^ I i i i i I i i i i I i i i i_ 

0.3 0.4 0.5 0.6 



n, (fm- 3 ) 



FIG. 2: Dimensionless parameters ^urii/^ (i = n ,A) versus n& in superfluid nucleon-hyperon 
matter. 

are different, we show the dimensionless combinations Ci-j^-i/ €2 as functions of the baryon 
number density rib- 

What is the nature of the relations (!83l) - (l85|) ? The relations ( |83|) follow from the equi- 
librium condition ([8]) with respect to the fast reaction ([7]). This condition is valid even if we 
allow for the dissipative processes dl])-®. Consequently, we can write 

Re(/i S ) + Re(/i p ) = Re(/i n ) + Re(/i A ) (86) 

or, in view of Eq. flTOj) (we remind that Re(/ij) = /ijdi S ; xt, p = xe + 

XEp = + ^A- (87) 

Then, substituting Eqs. (T4TT) (|49|) into Eq. (!H7|) and equating prefactors in front of the same 
four-divergences, we obtain the relations fl83l) . 

It is convenient to note that the constraint (1H7I) leads to a simple equation interrelat- 
ing four- vectors wfcy It holds true only if the bulk viscosities are generated by the non- 
equilibrium reactions. To derive this equation, we sum up the potentiality conditions 
(with proper signs) for all baryon species. Then, using Eqs. ([8]) and (IH7I) . we get 

= & [w\ n) + w[ K) - ufo - wfa] . (88) 



d" 



w 



(n) 



+ W 



(A) 



— W 



w 



(p) 
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In view of this equation and Eq. (T4"5l) , there are only 4 — 2 = 2 independent four- vectors 
in the system. 

Now let us explain why the bulk viscosities satisfy Eqs. (I8~4"j) and (I8"5j) . For that purpose, 
we consider the entropy generation equation. Neglecting all the dissipative processes except 
for the non-equilibrium reactions (jTJ) — (e.g., neglecting the thermal conductivity and shear 
viscosity), we can obtain from the hydrodynamics discussed in this section, 

Td^S" = 5fiAT = X5fi 2 . (89) 

A similar expression is valid for npe matter (see Eq. (79) of Ref. [l|). The chemical potential 
disbalance Sfj, is given by Eq. (1591 . Substituting into Eq. (I59"|) variations <Sn&, 5nn, and Sn-sn 
calculated from the continuity equations (I56 i) -( l58l) . we find 

1 \d5a ... 35a d5fi ,. . _ . 

fy = x /q fl r /q \ — ^— div(J 6 ) + - — div{J H ) + - div(J Sn ) . 

A {85fi/8n H — 85fi/8n-z n ) — iu> \_8n b 8n H 8n^ n 

(90) 

It follows from Eq. (|90|) . that for any given it is always possible to choose four- vectors 
in such a way that 5[i — at some point and in some particular moment (even if 
some baryon species are non-superfluid). In other words, this means that we can vanish the 
entropy generation rate (1891 at this point. 

On the other hand, the entropy generation equation in terms of the effective bulk viscosi- 
ties takes the form (see Eqs. (JHSJ and (I4"0l-(l4l 



[^e q df 1 [Y qk w^ k) ) + Cmpd^f d^Y^w^). (91) 



The right-hand side of this equation satisfies two conditions. First, it is a positive-definite 
quadratic form (entropy cannot decrease!). Second, according to Eqs. (|89|) and (J9~D]) . one can 
vanish it by an appropriate choice of w' 4 and (at some particular moment and at some 
point). There is a mathematical theorem that these two conditions are consistent with each 
other if and only if the determinant of the matrix composed of coefficients of the quadratic 
form vanishes. This result is independent of an actual number of superfluid baryon species. 
That is, the determinant will be zero in the case when all four baryon species are superfluid 
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as well as in the case when some of them are normal (for example, nucleons are normal, 
hyperons are superfluid, or neutrons are superfluid, while other particles are normal). As a 
consequence, we arrive at the six additional constraints (1511) and (1551) on the bulk viscosity 
coefficients. 

In this section we have assumed that baryons of all species are superfluid. However, the 
hydrodynamics formulated here can be easily extended to the situation when some baryon 
species are normal. In this case, one should vanish matrix elements related to these 
baryon species. For example, if neutrons are non-superfluid, then Y nk = Y kn = 0. 



IV. SUMMARY 



In this paper we have analyzed the bulk viscosity due to non-equilibrium particle trans- 
formations in superfluid nucleon-hyperon matter of neutron stars. Our approach is similar 
to that used in Ref. [lj for the case of superfluid npe matter. 

Our main results are as follows: 

(i) We have demonstrated that the expression for the bulk viscosity of normal (non- 
superfluid) nucleon-hyperon matter, widely used in the literature, is inaccurate. We have 
presented the correct derivation of the bulk viscosity and compared it with the results of Ref. 



151 ] . Numerically, both formulae give almost similar results (within a few tens of percent). 

(it) We have extended the hydrodynamics of superfluid mixtures reported in Refs. [l|, Q] 
to allow for a possible presence of superfluid hyperons. We have determined the general form 
of dissipative terms entering the equations of this hydrodynamics and showed that generally 
(when all baryon species are superfluid), it contains sixteen bulk viscosity coefficients. 

(Hi) We have calculated and analyzed the sixteen bulk viscosity coefficients assuming 
they are generated by non-equilibrium reactions (!Tl)-(j4"]) of particle mutual transformations. 
We have shown that only three of them are independent. All other 13 bulk viscosities 
can be expressed through these three using Eqs. (1501) and (l83l - (l85l) . In addition, we have 
demonstrated that Eq. (1671) for the bulk viscosity coefficient £2 formally coincides with Eq. 
(122]) for that in normal matter (however, the reaction rates (/ = 1, . . . , 4) are affected by 
superfluidity). 

Our results can be important for the studies of dynamical instabilities in pulsating su- 
perfluid neutron stars, especially for the studies of the r-mode instability. They can also be 
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important for modeling of the thermal evolution of pulsating neutron stars and for analyzing 
rotochemical and gravitochemical heating of millisecond pulsars with superfluid nucleon- 



hyperon cores 



in Refs. 
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for non-superffuid npe matter of neutron stars, this problem was considered 
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Let us notice, that to start such an analysis one needs to know the matrix Y^, which is 
the most important ingredient in the hydrodynamics of superfluid nucleon-hyperon mixture. 
To our best knowledge, it has not been considered in the literature so far. We are planning 
to fill this gap and present its extensive calculations in a subsequent publication. 
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